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ABSTRACT : 

In the last chapter of his book " The Algebraic Theory of Modular Systems " published in 1916, 
F. S. Macaulay developped specific techniques for dealing with " unmixed polynomial ideals " by 
introducing what he called " inverse systems " . The purpose of this paper is to extend such a 
point of view to differential modules defined by linear multidimensional systems, that is by linear 
systems of ordinary differential (OD) or partial differential (PD) equations of any order, with any 
number of independent variables, any number of unknowns and even with variable coefficients. 
The first and main idea is to replace unmixed polynomial ideals by pure differential modules. 
The second idea is to notice that a module is 0-pure if and only if it is torsion-free and thus if and 
only if it admits an " absolute parametrization " by means of arbitrary potential like functions, or, 
equivalently, if it can be embedded into a free module by means of an " absolute localization " . 
The third idea is to refer to a difficult theorem of algebraic analysis saying that an r-pure module 
can be embedded into a module of projective dimension r, that is a module admitting a projective 
resolution with exactly r operators. 

The fourth and final idea is to establish a link between the use of extension modules for such a 
purpose and specific formal properties of the underlying multidimensional system through the use 
of involution and a " relative localization " leading to a " relative parametrization " . 
The paper is written in a rather effective self-contained way and we provide many explicit examples 
that should become test examples for a future use of computer algebra. 
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1) INTRODUCTION : 



Let D = K[di, dn] — K[<]\ be the ring of differential operators with coefficients in a differen- 
tial field K with n commuting derivations di, 9„ and commutation relations diO — adi+dia, Va € 
K. If y^, ...,y™ are m differential indeterminates, we may identify Dy^ -\- ... -I- Dy™" — Dy with 

and consider the finitely presented left differential module M with presentation D™ 
M — > determined by a given linear multidimensional system with n independent variables, 
m unknowns and p equations. Applying the functor hom£i{», D), we get the exact sequence 
—i' homo{M, D) — D^' — > N — > of right differential modules that can be transformed 

by a side-changing functor to an exact sequence of finitely generated left differential modules. 
This new presentation corresponds to the formal adjoint ad{V) of the linear differential operator 
T> determined by the initial presentation but now with p unknowns and m equations, obtaining 
therefore a new finitely generated left differential module N and we may consider homi:i{M, D) 
as the module of equations of the compatibility conditions (CC) of ad{T>), a result not evident 
at all (Compare to [24]). Using now a maximum free submodule — > — > homo{M,D) 
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and repeating this standard procedure while using the well known fact that ad{ad{T>)) = V, 
we obtain therefore an embedding homD{hoinD{M, D), D) of left differential mod- 

ules for a certain integer \ < I < m because X is a field and thus D is a noethcrian bi- 
module over itself, a result leading to / = rkD{homD{M, D)) = rkoiM) < m as in ([19], p 
178,201)(See section 3 for the definition of the differential rankrko)- Now, the kernel of the map 
e : M — ;> homD{hom£){M,D),D) : m e{m){f) = f{m),'^f € homoiM^D) is the torsion sub- 
module t{M) C M and e is injective if and only if M is torsion-free, that is t{M) = 0. In that case, 
we obtain by composition an embedding — >• M — >^ D' of M into a free module that can also be 
obtained by localization if we introduce the ring of fractions S~^D = DS~^ when S = D — {0}. 
This result is quite important for applications as it provides a (minimum) parametrization of the 
linear differential operator V and amounts to the controllability of a classical control system when 
n = 1 ([24], p 258). This parametrization will be called an ''absolute parametrization " as it only 
involves arbitrary "potential-like " functions (See [1], [9], [18], [19], [20], [24], [25] and [32] for more 
details and examples, in particular that of Einstein equations). 

The purpose of this paper is to extend suh a result to a much more general situation, that 
is when M is not torsion-free, by using unexpected results first found by F.S. Macaulay in 1916 
through his study of " inverse systems " for " unmixed polynomial ideals " . 

For this we define the purity filtration : 

= tn{M) C tn-liM) C ... C ti{M) C toiM) = t{M) C M 

by introducing tr{M) = {m £ M \ cd{Dm) > r} where the codimension of Dm is n minus 
the dimension of the characteristic variety determined by m in the corresponding system for one 
unknown. The module M is said to be r-pure if tr{M) = 0,tr-i{M) = M or, equivalently, if 
cd{M) = cd{N) = r,yN c M and a torsion-free module is a 0-pure module. Moreover, when 
K = k = cst{K) is a field of constants and m = 1, a pure module is unmixed in the sense of 
Macaulay, that is defined by an ideal having an equidimensional primary decomposition. 

Example 1.1 : As an elementary example with K = k = Q,m = l,n — 2,p — 2, the differential 
module defined by d22y = 0, di2y = is not pure because z' = d2y satisfies d2z' = 0, diz' = while 
z" = diy only satisfies ^22" = and {{x2)'^ , X1X2) = (xi) ^ (Xi)X2)^- We obtain therefore the 
purity filtration = <2(M) C ti(M) C toiM) = t{M) = M with strict inclusions as 7^ 2;' e ti{M) 
while z'' e toiM) but z" ^ ti{M). 

Prom the few (difficult) references ([1], [9], [15], [18], [27]) deahng with extension modules ext^{M) = 
ext^j-,{M, D) and purity in the framework of algebraic analysis, it is known that M is r-pure if and 
only if there is an embedding — >^ M — >■ ext^jj{ext^jj{M, D), D). Indeed, the case r = is exactly 
the one already considered because ext%{M, D) = homoiM, D) and the ker/coker exact sequence: 

— > ext\N) — > M — > ext°{ext°{M)) — > ext^{N) — > 

allows to test the torsion-free property of M in actual practice by using the double-duality formula 
t{M) = ext^{N) as in ([19]). Also, when r > 1, a similar construction that we shall recall and 
illustrate in section 4 provides a finitely generated module L with projective dim,ension pduiL) = r, 
that is a minimum resolution of L with only r operators, and an embedding ^ M — )• L that 
allows to exhibit a relative parametrization of V because now the parametrizing potential-like func- 
tions are no longer arbitrary but must only depend on arbitrary functions of n — r variables. 

Example 1.2 : With K = k = Q,m = 2,n = 3,r = 1, the differential module M defined by 
the involutivc system <I>^ = d^y^ = 0, $^ = d^y"^ = 0, = d2y^ — diy"^ = is 1-purc and admits 
the resolution — > D — > — > — > M — > 0. The differential module L defined by 
the system d^z = is also 1-pure and admits the resolution — > D — > D — > L — > 0. We 
finally obtain the relative parametrization y^ = d\z, y'^ = d2Z providing the strict inclusion M C L. 

In a simple way, this result can be considered as a measure of how far a module is from being 
projective, recalling that a module P is projective if there exists another (projective) module Q 
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and a free module F such that P Q) Q F. 

We adapt the' " relative localization" technique used by Macaulay and combine it with the " in- 
volution" technique used in the formal theory of systems of partial differential equations in order 
to obtain an explicit procedure for determining L when M is given. Many examples will illus- 
trate these new methods that avoid the previous abstract arguments based on " double duality" . 
In particular, original non-commutative examples will also be presented. However, we point out 
the fact that the latter method can be adapted without any change to the case of systems with 
variables coefficients as it only depends on the use of adjoint operators but the following example 
will explain by itself the type of difficulty involved. 

Example 1.3 : Starting now with K = Q(x^, x^), m — 2,n — 3,r — 1, the new differential module 
M defined by d^y^ = 0, d^y^ = 0, d2y^ — diy"^ + x^y'^ = is also 1-pure and the differential module 
L is again defined by d^z = as in the previous example. However we obtain the totally different 
relative parametrization y^ = d\2Z — x'^d2Z + z,y'^ = ^22^ providing the strict inclusion M C L. 
More generally, we may consider a constant parameter a G k = Q and consider the new system 
day^ = 0, d^y"^ = 0, d2y^ — diy^ + ax^y"^ = depending on a. For a = we find back the case of 
the previous example and we let the reader wonder why the situation only changes when a 7^ 0. 

The content of the paper is just following the introduction. 

In section 2 we recall the definitions and results from the formal theory of systems of OD/PD 
equations that will be crucially used in the sequel. We pay a particular emphasize to the definition 
of involution and the way to introduce the Spencer operator in this framework. We also study the 
possibility and difficulty to use computer algebra in this framework. 

In section 3 wc recall the basic tools needed from module theory and homological algebra in 
a way adapted to our purpose, in particular the definition of the extension modules, and provide 
a few of their properties which, though well known by specialists of algebraic analysis, cannot be 
found easily in the literatiu-c. Meanwhile, we provide a few links with the preceding section which 
are not so well known. Many explicit examples will illustrate the main concepts in the commutative 
(constant field k) and the non-commutative (differential field K) framework. 

In section 4 we shall recall the proof of the theorem already quoted showing how to embed an 
r-pure module M into another module L with projrective dimension equal to r. We shall provide 
for the first time explicit computations of this result in order to point out the difficulty encountered 
in such a procedure as a motivation for avoiding it. 

In section 5 we extend the work of Macaulay, showing why only pure modules can fit with 
relative localization in a coherent way with what happens for torsion-free modules. Meanwhile, we 
shall extend for the first time this work to the non-commutative framework, showing in particular 
that the operator introduced by Macaulay ([11], §60) for studying inverse systems is nothing else 
than the Spencer operator. Many explicit examples, including highly non-trivial ones provided by 
Macaulay himself, will be fully treated in such a way that any engineer, even with a poor knowl- 
edge of homological algebra, will nevertheless become intuitively able to understand and apply 
these new techniques without reading the previous sections, just comparing to the way the same 
examples have been treated in section 4 by means of another approach. 

2) TOOLS FROM SYSTEM THEORY : 

If X is a manifold of dimension n with local coordinates (.t) = (.x^, ....t"), we denote as usual by 
T = T{X) the tangent bundle of X, by T* = T*{X) the cotangent bundle, by A''T* the bundle of 
r-forms and by SqT* the bundle of q-symmetric tensors. More generally, let S be a vector bundle 
over X , that is (roughly) a manifold with local coordinates (.t% y*) for i ~ 1, ...,n and fc = 1, m 
simply denoted by {x,y), projection it : E ^ X : (x,y) ^ (x) and changes of local coordinates 
X = (p{x), y = A{x)y. If E and F are two vector bundles over X with respective local coordinates 
{x,y) and {x,z), we denote by ExxF the fibered product of E and F over X as the new vector 
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bundle over X with local coordinates {x,y,z). We denote hy f : X E : (x) {x,y = f{x)) a 
global section of E, that is a map such that tt o / = zdx but local sections over an open set U <Z X 
may also be considered when needed. Under a change of coordinates, a section transforms like 
f{(f{x)) = A{x)f{x) and the derivatives transform like: 

^(^ix))d.cp-ix) = {d,Ai{x))f''ix)+A[ix)dj'{x) 
We may introduce new coordinates (x%y'^,2/*^) transforming like: 

We shall denote by Jq{E) the q-jet bundle of E with local coordinates (x', y*^, , y^^^, ...) = {x,yq) 
called jef coordinates and sections fq : (x) (x, /'^(x), /f (x), ...) = (x, /^(x)) transforming 
like the sections jq{f) : (x) ^ (x, /''(x), 5i/''(x), 9^ /''(x), ...) = {x,jq{f){x)) where both /g and 
are over the section / of E. Of course Jq{E) is a vector bundle over X with projection tt^ 
while Jq+r{E) is a vector bundle over Jq{E) with projection 7r|"'"'',Vr > 0. 

DEFINITION 2.1: A linear system of order g on is a vector sub-bundle Rq C Jq{E) and a 
solution of i?g is a section f of E such that is a section of 

Let fj, = (/xi, ...,fj,n) be a multi-index with length = /ii + ... + fJ-m class i if = ... = = 
7^ and /x + = {in, m-i, in + ...,/x„). We set = < /c < w, < < g} 

with = y'^ when =0. If S is a vector bundle over X with local coordinates (x*,?/*^) for 
i = l,...,n and fc = l,...,m, we denote by Jq{E) the g-jet bundle of with local coordinates 
simply denoted by {x,yg) and sections fq : (x) — > (x, /'^(x), /|^(x), /j^ (x), ...) transforming like 
the section : (x) — )■ (x, /'^(x), 9j/^(x), 9ij/^(x), ...) when / is an arbitrary section of E. 

Then both fq G Jq{E) and G Jq{E) arc over f £ E and the Spencer operator just allows 

to distinguish them by introducing a kind of "difference' through the operator D : Jq^i{E) 
T* ® Jq{E) : ^ - with local components (9i/^(x) - /f (x), 9i/j=(x) - 4(x),...) 

and more generally j(x) = dif^{x) — f^+i.{x). In a symbolic way, w/ien changes of 

coordinates are not involved, it is sometimes useful to write down the components of D in the 
form di = di — 5i and the restriction of D to the kernel Sq+iT* (g) £^ of the canonical projection 
7^^+^ : Jq+i(E) Jq{E) is minus the Spencer map 5 = dx' A (5, : S,j+iT* ® E ^ T* ® S'gT* ® E. 
The kernel of D is made by sections such that /g+i = ji{fq) = j2{fq-i) = ■■■ = jq+iif)- Fi- 
nally, if Rq C Jq{E) is a system of order q on E locally defined by linear equations ^'^{x,yq) = 
al.^{x)y^ =- and local coordinates (x, z) for the parametric jets up to order q, the r- prolongation 
Rq+r = Pr{Rq) = Jr{Rq) H Jq+r{E) C Jr{Jq{E)) IS locally defined when r = 1 by the lin- 
ear equations $'^(x,yg) = 0, rfi$'^(x, y^+i) = 0^^(2^)2/^+1. + 9ial.'^{x)y^ = and has symbol 
gq+r = Rq+r H Sq+r-T* ® E d Jq+r{E) if onc looks at thc top order terms. If G Rq+i is 
over fq ^ Rq, differentiating the identity al.^{x)f^{x) = with respect to x* and substracting the 
identity ii]f'{x)f'^j^i. {x)+diaY'{x)f'^{x) = 0, we obtain the identity a].^ {x){dif'^{x) — f'^j^i.{x)) = 
and thus the restriction D : Rg+i ^ T* i?, ([17], [18], [30]). 

DEFINITION 2.2: Rg is said to be formally integrable when the restriction ir'^Xr^^ ' ^g+r+i 
Rg+r is an epimorphism Vr > or, equivalently, when all the equations of order q + r are obtained 
by r prolongations only Vr > 0. In that case, Rg+\ C Ji{Rq) is a canonical equivalent formally 
integrable first order system on Rg with no zero order equations, called the Spencer form. 

Finding an intrinsic test has been achieved by D.C. Spencer in 1970 ([30]) along coordinate 
dependent lines sketched by M. Janet in 1920 ([7]) and W. Grobner in 1940 ([4], [6]). The key 
ingredient, missing in the old approach, is provided by the following definition. 

Let T* be the cotangent vector bundle of 1-forms on X and A^T* be the vector bundle of s-forms 
on X with usual bases {dx^ = dx^^ A ... A dx^^ } where we have set I = {ii < ... < is). Moreover, in- 
troducing the exterior derivative d : A^T* — > j-* . ^ _ uj{x)dx^ — > du = diU)i{x)dx^ A dx^ , 
we have = dod = and may introduce the Poincare sequence: 
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A°r* A A^T* A A^r* A ... A a"t* — ^ o 

PROPOSITION 2.3: There exists a map S : A^T* ® Sg+iT* (g) E ^ a^+^T* O SgT* (g) which 
restricts to 6 : A^T* (g) g^+i a'+^T* (g) and 5^ ^ 5 o (5 = 0. 

Proof. Let us introduce the family of s-forms co = {w^ = v^ ^dx^} and set {5u))^ = dx^ Aco^_^_i.. 
We obtain at once {S'^lo)^ = dx^ A dx^ A w^+i . = 0. 

Q.E.D. 

The kernel of each S in the first case is equal to the image of the preceding 6 but this may no 
longer be true in the restricted case and we set: 

DEFINITION 2.4: We denote by H'^+r{9q) cohomology at A'T* ® g^+r of the restricted 
5-sequence which only depends on Qq. The symbol Qq is said to be s- acyclic if -ffg+r = •.• = ^fg+r = 
0,Vr > 0, involutive if it is n-acyclic and finite type if gg+r = becomes trivially involutive for r 
large enough. 

DEFINITION 2.5: Rq is said to be involutive when it is formally integrable and its symbol gq 
is involutive, that is to say all the sequences ... -4 A^T* ® gg+r ••• are exact VO < s < n, Vr > 0. 

Equivalently, the following procedure, where one may have to change linearly the independent 
variables if necessary, is the heart towards the next effective definition of involution. It is intrinsic 
even though it must be checked in a particular coordinate system called 6-regular ([17], [18], [29]) 
and is particularly simple for first order systems without zero order equations. 

• Equations of class n: Solve the maximum number of equations with respect to the jets of 
order q and class n. Then call {x^, ...,x") multiplicative variables. 



• Equations of class i > 1: Solve the maximum number /J* of remaining equations with respect 
to the jets of order q and class i. Then call (a;-^, a;') multiplicative variables and ...,a;") 
non-multiplicative variables. 



• Remaining equations equations of order < q — 1: Call (a;^, ...,a;") non-multiplicative variables. 

In actual practice, we shall use a multiplicative board where the multiplicative "variables" are rep- 
resented by their index in upper left position while the non-multiplicative variables are represented 
by dots in lower right position. 

DEFINITION 2.6: A system of PD equations is said to be involutive if its first prolongation 
can be achieved by prolonging its equations only with respect to the corresponding multiplicative 
variables. In that case, we may introduce the characters a* = ^-^r^[^r^j^z^ ~ Pg for z = 1, n and 
we have dim{gg+\) = -|- ... -|- a". Moreover, one can exhibit the Hilbert polynomial dim{Rg+r) 
in r with leading term {a/(B?)r'^ with d < n when a is the smallest non-zero character in the case 
of an involutive symbol. Such a prolongation allows to compute in a unique way the principal 
{pri) jets from the parametric {par) other ones. This definition may also be applied to nonlinear 
systems as well. 

REMARK 2.7: For an involutive system with ^ = /3q < m, then {y^~^^ , ...,y"^) can be given 
arbitrarily and may constitute the input variables in control theory, though it is not necessary 
to make such a choice. In this case, the intrinsic number a = aq=m — (3>0 is called the 
n- character and is the system counterpart of the so-called "differential transcendence degree" in 

difFcrcntial algebra. As we shall see in the next section, the smallest non-zero character and the 
number of zero characters are intrinsic numbers that cannot be known without bringing the system 
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to involution and we have > ... > > 0. 

EXAMPLE 2.8: ([11], §38, p 40 where one can find the first intuition of formal integrabil- 
ity) The primary ideal q = ((xi)^)XiX3 ~ X2) provides the system yn = 0,yi3 —1/2 = which 
is neither formally integrable nor involutive. Indeed, we get d^yu — di{yi3 — 1/2) = 2/12 and 
dzy 12 — d2 (2/13 — 2/2) — V22i that is to say each first and second prolongation does bring a new second 
order PD equation. Considering the new system 2/22 = 0, 1/12 = 0, 1/13 — 2/2 = 0, i/n = 0, the question 
is to decide whether this system is involutive or not. One could use Janet or Grobner algorithm but 
with no insight towards involution. In such a simple situation, as there is no PD equation of class 3, 
two evident permutations of coordinates (1, 2, 3) (3, 2, 1) or (1, 2, 3) (2, 3, 1) both provide one 
equation of class 3, 2 equations of class 2 and 1 equation of clas 1. It is then easy to check directly 
that the first permutation brings the involutive system y-x^ = 0, 2/23 = 0, 1/22 = 0, yis — 2/2=0 that 
will be used in the sequel and we have al = 0, = 0, 0:2 = 2. 

EXAMPLE 2.9: With n = 4,m = l,q = 1, K = Q{x^,x'^,x^,x^), let us consider the system Ri: 

{2/4 - x^y2 - y = o, 2/3 - x'^yi = 

Again, the reader will check easily that the subsystem R[ C .Ri : 

u= 2/4 - x^yi - 2/ = 
v= 2/3 - x^Vi = 
w= 2/2-2/1=0 

namely the projection r\^^ of i?2 to Ri , is formally integrable and even involutive with one equa- 
tion of class 4, one equation of class 3 and one equation of class 2. 

In the situation of the last remark, the following theorem will generalizing for PD control sys- 
tems the well known first order Kalman form of OD control systems where the derivatives of the 
input do not appear ([27], VI, 1.14, p 802). For this, we just need to modify the Spencer form and 
we provide the procedure that must be followed in the case of a first order involutive system with 
no zero order equation, for example an involutive Spencer form. 

• Look at the equations of class n solved with respect to yj^,...,y^. 

• Use integrations by part like: 

2/^ - a{x)y^+' = dn{y' - a{x)y^+') + dr^a{x)y^+' =yi + dna{x)y^+^ 

• Modify y^,...,y^ to y^,...,y^ in order to "absorb" the various 1/^+^, ...,2/™ only appearing in the 
equations of class n. 

We have the following unexpected result providing what we shall call reduced Spencer form: 

THEOREM 2.10: The new equations of class n only contain i/f+\ with < i < n - 1 

while the equations of class 1, n — 1 no longer contain y^~^^, 2/™ and their jets. Accordingly, 
as we shall see in the next section, any torsion element, if it exists, only depends on y^, y^. 

Proof. The first assertion comes from the absorption procedure. Now, if 2/™ or 2/™ should ap- 
pear in an equation of class < n—l, prolonging this equation with respect to the non- multiplicative 
variable x" should bring or and (here involution is essential) we should get a linear com- 
bination of equations of various classes prolonged with respect to x^,...,x"'~^ only, but this is 
impossible and we get the desired reduced form. 

Q.E.D. 

When Rq is involutive, the linear differential operator V : E ^ Jq{E) Jq{E)/Rg — Fq of order 
q with space of solutions O C .B is said to be involutive and one has the canonical linear Janet 
sequence ([17], p 144): 

^ e ^ T A Fo A Fi A ... ^ ^ 
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where each other operator is first order involutive and generates the compatibility conditions (CC) 
of tfie preceding one. As the Janet sequence can be cut at any place, the numbering of the Janet 
bundles has nothing to do with that of the Poincare sequence for tlie exterior derivative, contrary 
to what many physicists beUeve. Moreover, the dimensions of the Janet bundles can be computed 
at once inductively from the board of multiplicative and non-multiplicative variables that can be 
exhibited for V by working out the board for Vi and so on. For this, the number of rows of this 
new board is the number of dots appearing in the initial board while the number nb{i) of dots in 
the column i just indicates the number of CC of class i for i = 1, n with nb{i) < nb{j),\/i < j. 
It follows that the successive first order operators Vi, ...,X>„ are automatically in reduced Spencer 
form. 

EXAMPLE 2.11: Coming back to Example 2.9 and changing slightly our usual notations, we 
get for X>i the following first order involutive system of CC in reduced Spencer form: 



(p^ = diV - d^u + x*diu - x^diV -v = 12 3 4 
(p"^ = d^w — d2U + diu — x^diw — w — Q 12 34 
(f?- = d^w — d2V + div — x'^diw = 1 2 3 • 

as d4U does not appear in and while u does not appear in 0^. 
We finally obtain for 1)2 the only CC: 

Tp = di<j)^ - d3<j)'^ - di(j)^ + x'^dicj)^ - x^di(j3^ - 0^ = 

DEFINITION 2.12: The Janet sequence is said to be locally exact at Fr if any local section of 

Fr killed by T>r+i is the image by Vr of a local section of ^V-i- It is called locally exact if it is 
locally exact at each F^ im < r < n. The Poincare sequence is locally exact, that is a closed 
form is locally an exact form but counterexamples may exist ([18], p 373). 

Equivalently, we have the involutive first Spencer operator Di : Cq = Rq ^ Ji{Rq) — >■ 
Ji{Rq) I Rq+i ~ T* (g) Rq/5{gq+i) = Ci of Order one induced by D : Rq+i — )■ T* (g) Rg. Introducing 
the Spencer bundles Cr ~ AT* (g) Rq/S{A^~^T* ® gq+i), the first order involutive (r + \)-Spencer 
operator Dr+i : Cr Cr+i is induced by D : AT* (g) Rq+i A'^+^T* ® Rq : a ^ ^q+i 
da®^q + {—lya A D^q+i and we obtain the canonical linear Spencer sequence ([17], p 150): 

^ e ^ Co A Ci A C2 ^ ... ^ Cn 

as the canonical Janet sequence for the first order involutive system Rq+i C Ji{Rq). 

The canonical Janet sequence and the canonical Spencer sequence can be; connected by a 
commutative diagram where the Spencer sequence is induced by the locally exact central horizon- 
tal sequence which is at the same time the Janet sequence for jq and the Spencer sequence for 
Jq+i{E) c Ji{Jq{E)) ([17], p 153) but this result will not be used in this paper (See [5], [20], [22], [23] 
for more details on Cosserat and Maxwell equations, see ([16]-[21]) and in particular ([22], [23]) for 
applications to engineering and mathematical physics) . 

REMARK 2.13: We shall revisit Example 2.8 in order to explain the word canonical " that 
has been used in the previous definitions. For this, starting with the inhomogeneous system 
2/33 = w,2/i3 — J/2 = we obtain easily the following inhomogeneous involutive system with its 
corresponding board of multiplicative and non- multiplicative variables: 



$4 


= 2/33 


= u 


1 


2 


3 


$3 


= y23 


= diu — d^v 


1 


2 


• 


$2 


= 2/22 


= diiu — di^v — d2V 


1 


2 


• 


$1 


= 2/13 - 2/2 


= V 


1 


• 


• 



Using prolongation with respect to the 4 non-multiplicative variables involved should bring 4 first 
order CC for the right members and we could wait for 4 third order CC involving u and v. 
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Surprisingly, we need the only CC \I/ = dssv — di3U + d2U = and obtain the differential sequence: 

as a single CC has no CC for itself (See ([18],p365) for the effective general procedure). 
Such a differential sequence is quite different from the canonical Janet sequence: 

which is the only sequence that can provide the Spencer sequence as we already said and could not 
be obtained by simply using Grobner bases. This remark will become essential in mathematical 
physics (foundations of continuum mechanics, gauge theory, general relativity) where only involu- 
tive operators must be used ([20], [22], [23]). We also check that the Euler-Poincare characteristic, 
namely the alternate sum of the circled dimensions of the vector bundles involved, does not depend 
on the differential sequence used as we get 1 — 2 + 1 = 1 — 4 + 4—1 = (See [18], p 378). 
In the same spirit, using certain parametric jet variables as new unknowns, we may set = y,z'^ = 



equation: 



y2,z = t/3 in order to obtain the following involutive first order system with no zero order 



class 3 
class 2 
class 1 



dsz'^ - = 0, ^3^2 - = 0, d^z^ = 0, dsz^ = 
d2Z^ -z^ = 0, d2Z^ - diz^ = 0, d2Z^ = 0, d2Z^ 
diz^-z'^ = 0,diz'^-z^ = 




where we have separated the classes. Contrary to what could be believed, this operator does not 
describe the Spencer sequence that could be obtained from the previous Janet sequence. Indeed, in- 
troducing the trivial vector bundle E with local coordinates (x^ , x'^,x^ ,y), it follows that Ji{E) has 
local coordinates [x^ ,x'^ ,z^ ^z"^ ,z^ ,z'^). Now, the involutive system R2 C J2{E) C Ji(Ji(E)) 
with involutive symbol 52 C S2T* (g) E C T* (g) T* ^ E C T* (g) Ji{E) projects onto Ji{E) 
but dim{R2) = 6 because par{R2) = {y, yi, 2/2, ys, 2/12} while we have only 4 unknowns 
{z^ , z"^ , z^ , z"^) . Nevertheless, as R2 projects onto Ji{E), we may construct a canonical Janet 
sequence for this new system where the successive Janet bundles involved will be the Spencer 
bundles Cr = A^T* €5 Ji{E)/S{A'^~^T* ig) 172) with a shift by one step in the numbering of the 
bundles as now Co = Ji{E) and the successive operators are induced by the composition of the 
inclusion R2 C J2{E) with the Spencer operator D : J2{E) — > T* Ji{E) as in ([17],pl44,150) or 
([18],p356). In any case, it is essential to notice that, both in the canonical Spencer sequence and 
in the canonical Janet sequence, any intermediate operator can be constructed explicitely without 
knowing the previous ones. 

EXAMPLE 2.14: With m = l,n = 4, g = 2, one could treat similarly the involutive system: 
y44 = 0, y34 = 0, yaa = 0, y24 — yis = with one equation of class 4, two equations of class 3 and 
one equation of class 2. 

EXAMPLE 2.15: Coming back to the involutive system of Example 2.9 with variable coeffi- 
cients, we let the reader prove that the Janet sequence is: 







e^0A@A@A0^o 



EXAMPLE 2.16: Let us finally consider the following involutive system of PD equations with 
two independent variables {x^^x"^) and three unknowns {y^,y'^,y^), where again a is an arbitrary 
constant parameter and we have set for simplicity y^ — diy^: 





+ y? 


+ y2 - 


y\ 


- ay^ 


= 


1 


2 


y\ 


-vl 


-2/1- 


y\ 


-ay^ 


= 


1 


2 


y\ 


-vl 


-2y? 






= 


1 


• 



Then the corresponding Janet sequence is: 

o^e^@A@A0^o 
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Finally, setting = y^ — y^,y^ = y^ + y^, we obtain the new first order involutive system: 



yj-yf-ay^ =0 1 2 
yl-yl-ay^ =0 12 

y\ -yf = o I 1 » 

with two equations of class 2 and one equation of class 1 in which y^ surprisingly no longer appears. 

If Xit-'tXti arc n algebraic indeterminatcs or, in a more intrinsic way, if x = Xidx^ £ T* is 
a covector and 2? : £; — Y F : ^ — y al''{x)d^^'' (x) is a linear involutive operator of order q, 
we may introduce the characteristic matrix a{x,x) = {a].'^ {x)xtJ. I M = and the resulting map 
cr^iP) '■ E — > F is called the symbol of T) at %. Then there are two possibilities: 

• If max^rk{a^{V) < m > 0: the characteristic matrix fails to be injective for any covector. 

• If max^rk{cr^{T>) = m <J=> a" = 0: the characteristic matrix fails to be injective if and only if 
all the determinants of the to x m submatriccs vanish. However, one can prove that this algebraic 
ideal a G if [x] is not intrinsically defined and must be replaced by its radical rad{a) made by all 
polynomials having a power in a. This radical ideal is called the characteristic ideal of the operator. 



DEFINITION 2.17: For each x ^ X, the algebraic set defined by the characteristic ideal is 
called the characteristic set of X> at a; and V = Ux^xVx is called the characteristic set of V. 

One has the following important theorem ([18], [29]) that will play an important part later on: 

THEOREM 2.18; (Hilbert-Serre) The dimension d{V) of the characteristic set, that is the maxi- 
mum dimension of the irreducible components, is equal to the number of non-zero characters while 
the codimension cd{V) = n — d{V) is equal to the number of zero characters, that is to the number 
of "full " classes in the board of multiplicative variables of an involutive system. 

EXAMPLE 2.19: Coming back to Remark 2.12, we obtain o = ((xa)^, X2X3, (X2)^, XiXs) =^ 
rad{a) = (x2,X3) and thus cd{V) = 2. However, if we take only into account Example 2.8, we 
should only get the radical ideal (xs) and the wrong result cd{V) = 1. The reason for using the 
radical can be seen from the equivalent first order system that shoul provide b = ((xa)^) •••) with ho- 
mogeneous polynomials of degree 4 and thus b C a with a strict inclusion though rad{a) = rad{b). 
A similar situation can be obtained with Examples 1.1 and 2.9. 



3) TOOLS FROM MODULE THEORY : 



We may roughly say that, if a reader familiar with Grobner bases ([4], [6]) and computer algebra 
looks at the previous section, he will feel cmbarassed because he will believe that " intrinsicness 
is always competing with complexity " as can be seen from Examples 2.8 -|- 2.12. However, even 
if he admits that it may be useful to have intrinsic and thus canonical procedures, then looking 
at the existing literature on differential modules ([1],[9],12]), he will really feel to be on another 
planet as the main difficulty involved in the theory of differentia modules is to understand why and 
where formal integrability and involution become essential tools to apply quite before dealing with 
the homological background of " algebraic analysis " involving extension modules. This is the main 
reason for which the case of variable coefficients is rarely treated " by itself" always refering to Weyl 
algebras for examples and the main difficulty we found when writing ([18], in particular Chapter 
IV). The central concept, essential for applications but well hidden in the literature dealing with 
filtred modules ([14] ,p 383) and totally absent from the use of Grobner bases because it amounts 
to formal integrability by duality, is that of a " strict morphism " . Accordingly, the purpose of this 
section will be to explain why such a definition, which seems to be purely technical, will be so 
important for studying extension modules and purity. 

If P = a'' dp G D ~ K[d], the highest value of with 7^ is called the order of the operator 
P and the ring D with multiplication (P, Q) — > P o Q = PQ is filtred by the order q of the 
operators. We have the filtration C K = Dq c Di c ... C Dq c ... C Dao = D. Moreover, it 
is clear that D, as an algebra, is generated hy K = Dq and T = Di/Dq with Di = K ®T ii we 
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identify an clement ^ — i^di E T with the vector field ^ — £^^{x)di of differential geometry, but 
with ^' G K now. ft follows that D — dDd is a bimodule over itself, being at the same time a left 
£)-module by the composition P — > QP and a right _D-module by the composition P — > PQ. 
We define the adjoint functor ad : D — > D°p : P = a^d^ — > ad{P) = (-l)l^ld^a'' and we have 
ad{ad{P)) = P. It is easy to check that ad{PQ) = ad{Q)ad{P),yP, Q G D. Such a definition can 
also be extended to any matrix of operators by using the transposed matrix of adjoint operators 
(See [18], [19], [22] for more details and applications to control theory and mathematical physics). 

Accordingly, ii y = (j/^, j/™) are differential indeterminates, then D acts on y'^ by setting 
diU'' = Ui — > di_i,y'' = y^ with diy^ = Vpt+ii ^^'^ Vo =2/*^- We may therefore use he jet coordi- 
nates in a formal way as in the previous section. Therefore, if a system of OD/PD equations is 
written in the form = a^^y^ ~ with coefficients a G K, we may introduce the free differential 
module Dy = Dy^ + ... + Dy"^ ~ and consider the differential submodule / = -D$ C Dy 
which is usually called the module of equations, both with the differential module M = Dy/D^ 
or D-module and we may set M = nM if we want to specify the ring of differential operators. 
The work of Macaulay only covers the case m = 1 with K replaced by A; C cst{K). Again, we 
may introduce the formal prolongation with respect to di by setting di^"^ = a]j^y'^_^i. + diO^f'y'^ 
in order to induce maps di : M — > M : y^ — > y^+i. by residue if we use to denote the residue 
Dy — > M : — > y^ by a bar as in algebraic geometry. However, for simplicity, we shall not 
write down the bar when the background will indicate clearly if we are in Dy or in M . 

As a byproduct, the differential modules we shall consider will always be finitely generated 
{k = l,...,m < oo) and finitely presented (r = l,...,p < oo). Equivalently, introducing the 
matrix of operators V = {a]J^d^) with m columns and p rows, we may introduce the morphism 

DP ^ D-^ : (Pr) — > (fVS^) : P — > P^ = PV over D by acting with D on the left of these 
row vectors while acting with V on the right of these row vectors and the presentation of M is 
defined by the exact cokernel sequence D^ — > — > M — > 0. It is essential to notice that 
the presentation only depends on K, D and $ or 2?, that is to say never refers to the concept of 
(explicit or formal) solutions. It is at this moment that we have to take into account the results 
of the previous section in order to understant that certain presentations will be much better than 
others, in particular to establish a link with formal integrability and involution. 

It follows from its definition that M can be endowed with a quotient filtration obtained from 
that of D™ which is defined by the order of the jet coordinates yq in Dqy. We have therefore the 
inductive limit C Mo C Mi C ... C C ... C = M with djM, C M,+i and M = DMq for 
q^O with prolongations DrMq C Mq+r,'^q,r > 0. 

DEFINITION 3.1: ([14],p 383) If M and are two differential modules and / : M — ^ A'' is 
a morphism over D compatible with the filtration, that is if f{Mq) C Ng with induced morphism 
fq : Mq — > Nq, then / is a strict morphism if fq{Mq) = f{M) n Ng, Vq > 0. 

Equivalently, chasing in the following diagram: 



Mg 

i 

M 

then / is strict if the induced morphism coker{fg) — > coker{f) is a monomorphism \/q > 0. 

DEFINITION 3.2: An exact sequence of morphisms finishing at M is said to be a resolution of 
M. If the differential modules involved apart from M are free, we shall say that we have a free 
resolution of M. Moreover, a sequence of strict morphisms is called a strict sequence. 

LEMMA 3.3: If / is a strict morphism as in the last definition, there are exact sequences 





; i 

Nq — > coker{fg) 

1 i 

N — > coker( f) 
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— > coker{fq) — > cofcer(/,+i), Vg > 0. 



Proof. As / is compatible with the filtrations and Mq C Mq+i,Nq C Nq+i, we have an induced 
morphism coker{fq) — > coker{fq+i). Now, as / is also strict, we have the following commutative 
and exact diagram: 



— > coker{fq) — > coker{f) 

i II 
— !■ coker{fq^i) — coker{f) 

The lemma finally follows from an elementary chase. 



Q.E.D. 



Having in mind that K is a left D-module with the standard action {D, K) — > K : [di, a) — > 
diC and that D is a bimodule over itself, we have only two possible constructions: 

DEFINITION 3.4: We define the system R= homK{M, K) = M* and set Rg = homK{Mg, K) = 
M* as the system of order q. We have the projective limit R = R^o — > ••• — > Rq — > ••• — > 
Ri — > Rq. It follows that fq € Rq : — > f^ £ K with a^f.^ f^ ~ defines a section at order 
q and we may set /oo = f € R for a section of R. For an arbitrary differential field K, such a 
definition has nothing to do with the concept of a formal power series solution (care). 

DEFINITION 3.5: We may define the right differential module homD{M,D). 

PROPOSITION 3.6: When M is a left £>-module, then R is also a left D-module. 

Proof As D is generated by K and T as we already said, let us define: 

{af){m) = af{m), \fa€ K,\/me M 

In the operator sense, it is easy to check that dia = adi + dia and that ^rj — r]£_ — 77] is the 
standard bracket of vector fields. We finally get {dif)'j^ — {dif){y^) = 9;/^ — fj^+i^ and thus 
recover exactly the Spencer operator of the previous section though this is not evident at all. We 
also get {didjf)l = dijf^ - - Sj/^+i, + =^ didj = djdi,yi,j = 1, ...,n and thus 

diRq+\ C Rq diR C R induces a well defined operator R — > T* ^ R : f — > dx^ ® dif. This 
result has been discovered (up to sign) by Macaulay in 1916 ([11]). For more details on the Spencer 
operator and its applications, the reader may look at ([22], [23]). 

Q.E.D. 

As £) is a bimodule over itself, it follows from this proposition that that D* = homK{D, K) 
is a left D-module. Moreover, using Baer's criterion ([28]), it is known that D* is an injective 
D-module as there is a canonical isomorphism: 

M* = homK{M,K) ~ homD{M,D*) 
where both sides are well defined ([2], Prop 11, p 18;)([28], p 37). 

DEFINITION 3.7: With any differential module M we shall associate the graded module 
G = gr{M) over the polynomial ring gr{D) ~ K[x] by setting G = ©^o<^g '^it^ (^q = ^q/^q+i 
and we get gq = G* where the symbol gq is defined by the short exact sequences: 

Mq_l -^Mq^Gq^O =^ ^ gq ^ Rq ^ Rq-1 — )■ 

We have the short exact sequences — > Dq-i — > Dq — > SqT — > leading to grq{D) ~ SqT 
and we may set as usual T* = homK{T, K) in a coherent way with differential geometry. Moreover 
any compatible morphism / : M — )■ N induces a morphism gr{f) : gr{M) — > gr{N). 
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EXAMPLE 3.8: UK — Q(a;),TO ~ l,n = 1, let us consider the system jjxxx — IJx = for which 
we may exhibit the basis of sections {/ = (1, 0, 0, 0, ...), /' = (0, 1, 0, 1, ...), /" — (0, 0, 1, 0, ...)} as in 
([ll],§59,p 67) or ([21]). We obtain d^f = 0, d^f = -f-f^d^r = -f and check that all the sec- 
tions can be generated by a single one, namely /" which describes the power series of ch{x) — l. With 
now m = 2, let us consider the module defined by the system yl^ = 0, = 0. Setting y = y^ — xy^ , 

we successively get y^ = -xyl-y^^y^^x = -^yl,yxxx = t/^ = %yxx-yx,y'^ = \yxx-xyx+y 
and a differential isomorphism with the module defined by the new system yxxx = 0- All the sec- 
tions of the second system are easily seen to be generated by the single section / = (0.0, 1, ...), a 
result leading to the only generating section f^{x) = f , /i = —\, fix = 0) •••) /^(s^) = fx = 0, ... 
of the initial system but these sections do not describe solutions because dxf^ — /.^ = 1 7^ and 
dxf^ — fx = X ^ 0. We do not know any reference in computer algebra dealing with sections (See 
[21] for more details) 

Coming back to the presentation of M under study, we notice that the morphism V involved 
is not compatible unless we shift the index of the filtration by ordij)) = q. In that case, we obtain 

im{T>) = I C Dy and may set Iq+r = I r\ Dq-^-rV but we have in general D^^V C /g+r only, that 
is the equations of order q-\-r may not be produced by r prolongations only. We have thus obtained: 

PROPOSITION 3.9: The morphism induced by V is strict if and only if V is formally intc- 
grable. Accordingly, the module version of both the Janet sequence and the Spencer sequence are 
strictly exact sequences. 

Proof. Using q-\- r instead of q in Lemma 3.3 and applying homK{*,K), we obtain the epimor- 
phisms Rq+r+l — > Rq+r — > 0, Vr > 0. 

Q.E.D. 

The reader will find in ([18], IV, 3) more details on the relations existing between G and M 

which are needed in order to study the non-commutative situation, at least when is a dif- 
ferential field as such a case is hard enough. We obtain in particular the Hilbert polynomial 
dimK{Mqj^r) = dim K {Rq+r) = ^'''^ + ••• where a is called the multiplicity of M and we use to set 
cdoiM) = cd{M) =n-r, rkoiM) = rk{M) = a if cd{M) = and otherwise. 

EXAMPLE 3.10: Coming back to the Example 2.8 of Macaulay, we obtain from Remark 2.13 

the free resolution — > D — > — > D — > M — > but only the morphism on the left is strict 
as for the morphism on the right we know that its image is indeed I2 = {2/33, 2/23, 2/22, 2/13 — 2/2} and 
not just {2/33, 2/13 — 2/2}- However, bringing the system to involution, we get the strict free resolution 
— > D — > D'^ — > — > D — > M — > as the module version of the Janet sequence and we 
let the reader exhibit the module version of the corresponding Spencer sequence as an exercise. 

If M is a differential module over the ring D = K[d] of diflFcrential operators and m G M, then 
the differential submodule Dm C M is defined by a system of OD or PD equations for one unknown 
and we may look for its codimension cd{Dm). A similar comment can be done for any differential 
submodule M' c M. Sometimes, a single element m £ M, called differentially primitive element, 
may generate M if Dm = M. 

EXAMPLE 3.11: With K = Q, let us consider the differential moduleM defined by the system 
vlx -y^ =0,yl = 0- Then, setting y = y^-y'^,we get yx = yl, yxx = ylx = y^ ,yxx - y = y"^ and 
thus yxxx — 2/x = provides another way to describe M by means of a single element as in ([18], 
p435). We have the following commutative and exact diagram: 

— )■ — y £,2 — ^ ]^ — ^0 

i i II 

— D — > D — > M — >0 

(P,Q) {P{d^-l),Qd) 

{Pd + Q) {{Pd + Q){d-' -d)) 
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If now we consider the differential module M defined by j/^^ — ay^ = 0, = where a is a constant 
parameter, we cannot find a differentially primitive element when K = Qif a = but we can when 
K = Q{x) for any value of a, as in Example 3.8. 

We may check the following definition in a constructive way ([27]): 

DEFINITION 3.12: UiM) = {m € M | cd{Dm) > r} is the greatest differential submodule of 
M having codimension > r. 

PROPOSITION 3.13: cd{M) = cd{V) = r^ a^""" 7^ O.a^-'^+i = ... = = <^ tr{M) ^ 
M,tr-i{M) = ... = to{M) = t{M) = M and this intrinsic result can be most easily checked by 
using the Spencer form of the system defining M. 

We are now in a good position for defining and studying purity for differential modules. 

DEFINITION 3.14: M is r-pure 4=^ U{M) = 0,tr-i{M) = M <S=4> cd{Dm) = r,Vm € M. In 
particular, M is 0-pure if t{M) = and, if cd{M) = r but M is not r-pure, we may call M/tr{M) 
the pure part of M . It follows that tr-i{M)/tr{M) is equal to zero or is r-pure (See the picture in 
[18], p 545). Finally, when tr-i{M) = tr{M), we shall say that there is a gap in the purity filtration: 

= tn{M) C t„_i(M) C ... C ti{M) C to{M) = t{M) C M 
PROPOSITION 3.15: tr{M) does not depend on the presentation or on the filtration of M. 

EXAMPLE 3.16: If K = Q and M is defined by the involutive system ^33 = 0, 2/23 = 0, 2/13 = 0, 
then z — 2/3 satifies d^z = 0, d2Z = 0, diz = and cd{Dz) = 3 while z' = 7/2 only satisfies d^z' = 
and cd{Dz') = 1. We have the purity filtration = h\M) c t2{M) = ti{M) c to{M) = t{M) = M 
with one gap and two strict inclusions. 

We now recall the definition of the extension modules ext]^ (M, D) that we shall simply denote 
by exf{M) and the way to use their dimension or codimension. We point out once more that 
these numbers cannot be obtained without bringing the underlying systems to involution in order 
to get informations on M from informations on G. We divide the procedure into four steps that 
can be achieved by means of computer algebra ([27]); 

• Construct a free resolution of M, say: 

... — > Fi — > ... — > Fi — > Fq — > M — >0 

• Suppress M in order to obtain the deleted sequence: 

... ^Fi^...^Fi^Fo^O 

• Apply honiDi*, D) in order to obtain the dual sequence heading backwards: 

... i — homD(F;.D) < — ... < — homD{Fi,D) i — homD{Fo,D) i — 

• Define ext'^{M) to be the cohomology at homD{Fi,D) in the dual sequence with ext°{M) = 
homD{M,D). 

The following nested chain of difficult propositions and theorems can be obtained, even in the 
non- commutative case, by combining the use of extension modules and bidualizing complexes in the 
framework of algebraic analysis. The main difficiilty is to obtain first these results for the graded 
module G = gr{M) by using techniques from commutative algebra before extending them to the 
filtred module M as in ([1], [9], [18], [19]). 

THEOREM 3.17: The extension modules do not depend on the resolution of M used. 
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PROPOSITION 3.18: Applying homr){», D) provides right I?-niodules that can be transformed 
to left _D-modules by means of the side changing functor and vice-versa. Namely, if No is a right 
Z)-module, then dN ~ A"T0kNd is the converted left L'-module while, if dN is a left Z)-module, 
then Nd = /\'^T*®kdN is the converted right D-module. 

PROPOSITION 3.19: Instead of applying homD{», D) and the side changing functor in the 
module framework, we may use ad in the operator framework. Namely, to any operator V : E — > F 
we may associate the formal adjoint ad{V) : A"r* ® F* — > /\nrp* ^ £.* .^j^j^ ^j^g useful though 
striking relation rkoiadiV)) = rkoiT^)- 

PROPOSITION 3.20: ext'{M) is a torsion module VI < i < n but ext^{M) = homD{M,D) 
may not be a torsion module. 

EXAMPLE 3.21: When M is a torsion module, we have homD{M,D) = (exercise). When 
n = 3 and the torsion-free module M is defined by the formally surjective div operator, the formal 
adjoint of div is —grad which defines a torsion module. Also, when n = 1 as in classical control 
theory, a controllable system allows to define a torsion-free module M which is free in that case 
and homoiM, D) is thus also a free module. 



THEOREM 3.22: ext\M) = 0, Vi > n + 1. 

THEOREM 3.23: cd{ext'{M)) > i. 
PROPOSITION 3.24: ext'{M) = 0,Vi < cd{M). 

THEOREM 3.25: cd{M) >r<^ ext\M) = 0, Vi < r. 



PROPOSITION 3.26: cd{M) = r cd{exV{M)) = r and exf{M)is r-pure. 

PROPOSITION 3.27: exf" {exf" {M)) is equal to or is r-pure, VO < r < n. 



PROPOSITION 3.28: If we set t-i{M) = M, there are exact sequences VO < r < n: 

— > tr{M) — > tr-i{M) — > exe{exe{M)) 

THEOREM 3.29: If cd{M) = r, then M is r-pure if and only if there is a monomorphism 
— > M — y exV{exr{M)) of left differential modules. 

THEOREM 3.30: M is pure <^ ext' {ext' {M)) = 0,Vs ^ cd{M). 

The last two theorems are known to characterize purity but it is however evident that they are 
not very useful in actual practice. 

THEOREM 3.31: When M is r-pure, the characteristic ideal is thus unmixed, that is a finite 
intersection of prime ideals having the same codimension r and the characteristic set is equidimen- 
sional, that is the union of irreducible algebraic varieties having the same codimension r. 

REMARK 3.32: For the reader knowing more about commutative algebra, we add a few de- 
tails about the localization used in the primary decomposition of a module which are not so well 

known ([3], [18], [21], [31]). For simplicity, setting k = cst{K), we shall denote hy A = k[x] the 
polynomial ring isomorphic to D = k[d] and consider a module M over A. We denote as usual 
by spec{A) the set of proper prime ideals in A, by max{A) the subset of maximal ideals in A 
and by ass{M) = {p G spec{A)\30 m e A-I,p = anuAim)} the (finite) set {pi,...,p/,} of 
associated prime ideals, while we denote by {pi, ...ps} the subset of minimum associated prime 
ideals. It is well known that M ^ ass{M) ^ 0. We recall that an ideal q C A is p- 
primary if a6 G q, & ^ q ==> a G rad{c{) = p G spec{A). We say that a module Q is p-primary 
if am = 0,0 7^ rn G Q =4> o G p = rad(q) G spec{A) when q = annA{Q) or, equivalently, 
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ass{Q) = {p}. Similarly, we say that a module P is p-prime if am = 0,0 ^ rn, E P a € 
p G spec{A) when p = annA{P)- It follows that any p-prime or p-primary module is r-pure with 
n — r = trd{A/p), a result generalizing ([11], §4, p 43). Accordingly, a module M is r-pure if and 
only if a = annA{M) admits a primary decomposition a = qin...nqs and rad{a) = pin...nps with 
cd{A/pi) = cd{M) = r,\/i = 1, s. In that case, the monomorphism — > M — > ®peass(M)-^p 
induces a monomorphism — > M — > Qi® ■■■®Qs called primary em,hedding where the primary 
modules Qi are the images of the localization morphisms M — > Mp. = S~^M with S = A — p 
inducing epimorphisms M — > Qi — > for i = 1, ...,s. Macaulay was only considering the case 
M = A/a with primary decomposition a = qi fl ... fl q^. 

EXAMPLE 3.33: With A; = Q and n = 3, then a = rad(o) = (xi, X2X3) = (xi, X2) n (xi, Xs) is 
unmixed and M = A/ a is 2-pure while a = ra(i(a) = (xiX2,XiX3) = (Xi)l~l(X2, Xa) is mixed, though 
an intersection of two minimum prime ideals and M = A/a is not 1-pure. On the contrary, if one 
has the primary decomposition a = ((Xl)^ X1X2, XiX3, X2X3) (xi, X2) n (xi, X3) n (xi, X2, X3)^ = 
pi n p2 n and M = A/a, then ass{M) = {pi,p2,m} with p,; c m for i = 1,2, though 
rad{a) = pi n p2 as before. In that case, there is an embedding — > M — > Qi ® Q2 (B Q3 
where Qi = A/pi is the image of the localization morphism M — > Mp. for i = 1,2 because pi 
is killed by X3 G ^ — Pi) p2 is killed by X2 S ^ — p2 and Q3 = A/vx^ is m-primary because 
rad(m^) = m G max{A). We have also an embedding — > M — > Mp^ © Mp^ ® but no 
element of the multiplicative setA — m = {l-l-a|aem} can kill any element of M and the image 
of M into Mm is thus isomorphic to M which is not a primay module. It is important to notice 
that the example of Macaulay q = ((X3)^)X2X3) (X2)^,XiX3 ~ X2) provides a p-primary module 
A/q with p = (x3,X2) even though the annihilating ideal of G = gr{M) is the homogeneous 
ideal = ((X3)^, X2X3, (X2)^, X1X3) = ((X2)^,X3) n (xi,X2,X3)^ which is a mixed ideal because 
ass{A/a) = {(X2,X3), (xi.X2,X3)}- However, we get rad{a) = (X2,X3) in a coherent way. 

PROBLEM : Is it possible to have a test for checking whether a differential module is pure or 
not without using the previous results ? 

4) MOTIVATION : 

As we already said in the introduction and in the previous section, a torsion-free module M 
is 0-pure because in that case to{M) = t{M) = 0. Accordingly, M can be embedded into a free 
module F and the inclusion, which may not be strict when n > 1, provides a parametrization by 
means of a finite number of potential-like arbitrary functions in the classical language of elastic- 
ity (Airy function) or electromagnetism (EM 4-potential) . As it is clear that such a situation is 
only a very particular case of purity, it remains to wonder what can happen for an r-pure module 
whenever r > 1. One has the following result ([9], [18], compare to [1], p494): 

THEOREM 4.1: If M is an r-pure differential module with r > 1, there exists a differential 
module L with pd{L) < r and an embedding M C L. 

Proof. First of all we notice that we have r > and thus any element to G Af is surely a torsion ele- 
ment because cd{Dm) > 0, that is M = t{M) is a torsion module with ext'^{M) = homD{M, D) = 
because D is an integral domain. Let now ... — > Fr — > ... — > Fi — > Fq — > N — > be a free 
resolution of the right differential module N = ext^(M). According to Proposition 3.26, we have 
cd{N) = r > too and N is also a torsion module with ext^{N) = homD{N, D) = 0. Applying the 
functor hom,u{». D) to the previous sequence or , equivalently, constructing the adjoint sequence in 
the operator framework while using the fact that ext{N) = 0,Vi < r according to Theorem 3.25, 
we obtain the finite long exact sequence with exactly r morphisms because N is finitely presented 
and exV{N) ^ 0: 

— ^ homD{Fo,D) — ^ ... — > homD{Fr-i,D) — > h(miD{Fr,D) — >L — ^ 

where the left differential module L is the cokernel of the last morphism on the right. As 
homD{F,D) is free whenever F is free because of the bimodule structure oi D = dF>d, the 
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corresponding deleted complex is: 

— y homD{Fo, D) — > ... — > homD{Fr-i, D) — > homD{Fr, D) — > 

Applying again homD{», D) and using the reflexivity of any free module F, that is the isomorphism 
homi3{homD{F,D),D) ~ F, we obtain the dual sequence: 







F\ — >■ Fq 







and a similar procedure may be followed with operators as we shall see in the next illustrating 
examples ([18], [27]). This sequence is exact everywhere but at and at Fo where its cohomology 
is just N by definition, that is to say ext^{L) = N = ext^{M). Looking for the cohomology at 
homoiFr, D) in the sequence obtained by duality from the resolution of A'' with coboundry module 
Bj. and cocycle module Zj-, we obtain the following commutative and exact diagram: 



Finally, composing the bottom monomorphism with the monomorphism — > M — > exf{N) 
provided by Theorem 3.29, we get the desired embedding M C L. It must be noticed that such a 
procedure can be followed equally well in the commutative and non-commutative framework, that 
is when if is a field of constants or a true differential field. 

Q.E.D. 

EXAMPLE 4.2: With K = Q,m = l,n = 4,q = 2, let us study the 2-pure differential module 
M defined by the involutive system: 









; 


; 




= Bf 




; 




— y homoiF, 




i 


exV{N) 


— > L 


i 


i 









$4 


= 2/44 


= 


1 


2 


3 


4 


$3 


= 2/34 


= 


1 


2 


3 


• 


$2 


= 2/33 


= 


1 


2 


3 


• 


#1 


= 2/24 - 2/13 


= 


1 


2 


• 


• 



From the board of multiplicative variables wc may construct at once the Janet sequence: 

^ e 

where T>i is defined by the involutive system: 







^4 


= di^^ 




= 


1 


2 


3 


4 


^f3 


= di<P^ 




= 


1 


2 


3 


4 


^f2 


= d4<i>^ 


- d2<^^ + di'i'^ 


= 


1 


2 


3 


4 


^1 


= d3$^ 


-d2*^ + di$2 


= 


1 


2 


3 


• 



and I?2 by the (trivially) involutive system: 



We have therefore the resolution: 



12 3 4 







D 



D 



M 
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leading to pd{M) < 3 and the deleted complex is: 

— >D — >D'^ — yD'^ — >D 
Applying homD{*,D) to this sequence, we get the sequence: 

O^f — Di — D'^ i — D'^ i — D 
which can be described by the following adjoint sequence: 

which is not a Janet sequence. As M is a torsion module, using now Theorem 3.25 wc get 
ext^{M) = 0, ext^{M) = and we check that N = ext^{M) ^ 0. For this, dualizing * by A 
and Q. by 6, we have to look for the CC of the inhomogeneous system: 



r ^1 - 










= x^ 


■ ^3 _ 


di6 


= A3 


^4 _ 


-d2e 


= A4 



which are not already provided by the system: 



< 



r $1 

$3 
$4 



-d^X^ - dsX^ 
-diX^ - rfiA^ 

-(d4A4-rf2Ai) + (rf3A3-rfiA2) 

fisA^ + daA^ 










One can check that the torsion module N can be generated by {m = d2A^ + diA**, = dsA^ — diA^} 
satisfying the involutive system: 



4 


= diU — 


d\v 


= 


1 


2 


3 


4 




= diV 




= 


1 


2 


3 


4 


2 


= dsu - 


d2V 


= 


1 


2 


3 


• 


1 


= d^v 




= 


1 


2 


3 


• 



with the two CC: 



= did)^ 



^30"* + d24>^ - di^ 



,1 - 



= 



= di(j)^ - d3(j)^ = 

Accordingly, we have the following strict free resolution of N: 

— — — >D'^ — > N — > 

with deleted complex: 

— > — > — > — >0 
Applying homj:,{»,D), we get the desired resolution of L, namely: 



1 


2 


3 


4 


1 


2 


3 


4 















Dualizing ^ by 2;, we finally discover that L is defined by the involutive system: 



d4Z^ — diz^ 

diz"^ 

daz'^ - diz"^ 
dsz^ 



12 3 4 

12 3 4 

1 2 3 • 

12 3* 
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and is therefore 2-pure with pd{L) < 2 and a strict inclusion M G L defined by y = z^. 



REMARK 4.3: In this example, we discover that, if L were also r-pure, we should therefore 
have an embedding — > L — > exV {exf^ (L)) = ext^{N) and thus an isomorphism exV{N) = L 
leading to an isomorphism Zr = homniFr, D) and to ^V+i = 0, as can be checked on this example 
with r = 2. It has been a challenge for the author during many months to find the following 
counter-example showing that, sometimes L may not even he a torsion module. 

EXAMPLE 4.4: According to the proof of the theorem. A'' = exV{M) does not depend on the 
resolution of M used while L does indeed depend on the resolution of N used. Coming back to 
the system studied in Example 2.8 and Remark 2.12 with r = 2, wo may use the shortest finite 
free resolution of M already presented, namely — > D — > — > D — > M — y 0. Therefore, 
taking the adjoint of the only CC found, we may define A'' by the system: 



V 

-u 



dz^e 

di^e + d26 



and obtain the corresponding involutive system: 



4 


= d^^e 


= 


1 


2 


3 




= d23e 


= 


1 


2 


• 


2 


= d220 


= 


1 


2 


• 


1 


= di^e 


+ d2e =0 


1 


• 


• 



We obtain the first order involutive system of CC: 





= dz 


f~d2 




= 


1 


2 


3 




= ds 


f-d2 


f 


= 


1 


2 


3 




= ds 


f-di 






1 


2 


3 




= d2 


i>'-di 




^2 =0 


1 


2 


• 




d3^'- 


d2^2 + 


dii'" 


+ = 









with the only CC : 

Wc may therefore introduce in reverse order the corresponding adjoint operators of the ones in- 
volved in the Janet sequence we have just constructed: 



-di\ 

X 
d2\ 
-daX 



= 1/2 
- ,/l 



02 

-{dsi^ 
is 



—dsu^ — 



4 - + (d2l/3 



-d^v'^ - d2V^ 



= 
= 
= 
= 



This last operator is defining L but is not involutive. We have the two torsion elements: 

div^ + i/^ 



doU^ - ;y2 



which are generating ext^{N) and are easily seen to satisfy the involutive system: 



ds/j/' - IJ'' = 0, 



dsA^' = 0, 



d2tJ.^ - d^i^ = 0, 



d2ir' = 



because ^2/^^ = ^2/"'^ + c^sA*^ + dxij} = 0. Finally, using the first equation, we may eliminate iJ" 
and identify /^^ with y because we have indeed ^33/1^ = 0, di^fj,^ — d2H^ = in order to obtain the 
strict inclusion M C L. Equivalently, we may also eliminate and respectively from fi^ and 
fi^ in order to obtain: 



disidiiy"" + u^) - d2{diu'' + u^) = 0, 
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but wc may notice that L is not 2-purc and thus a torsion module because v'^ (similarly 7/^) is 
not by itself a torsion element of L. Such a situation is well known in control theory with the 
SISO (single input u, single output y) system y — ii = because u (similarly y) is not by itself a 
torsion element but z = y — u is a torsion element because i = (See the pages 9 and 10 of the 
introduction in [17] for more details on such a comment). 

PROBLEM : Is it possible to find an analogue of the previous theorem or of the case r = 0, 
where L should be also r-pure with a free resolution having exactly r morphisms ?. 

5) ABSOLUTE AND RELATIVE LOCALIZATIONS : 

Surprisingly, the; positive; answer to such a problem has been given by Macaulay in ([11]) for 
differential modules defined by systems with constant coefficients and only one unknown. Our pur- 
pose in this section is to generalize this resul to arbitrary differential modules defined by systems 
of PD equations with coefficients in a differential field. 

Now we hope that, after reading the previous section, the reader is convinced that the use of 
extension modules is a quite important though striking tool for studying linear multidimensional 
systems. Of course, as for any new language, it is necessary to apply it on many explicit examples 
before being familiar with it. However, it is evident that it should be even more important to have 
a direct approach allowing to exhibit the purity filtration and, in particular, to recognize whether 
a differential module is pure or not. The purpose of this section is to combine the module approach 
with the system approach, while taking into account the specific properties of the Spencer form in 
a way rather similar to the use of the Kalman form of a control system when testing controllabil- 
ity, namely to check that an ordinary differential module is 0-pure. For this, we shall divide the 
procedure into a few successive constructive steps that will be illustrated on explicit examples. 

• STEP 1: Whenever a system Rq C Jq{E) is given, there is no way to obtain informations on 
the corresponding module without bringing this system to an involutive or at least formally inte- 
grable system by means of prolongations and projections as in the Example 2.8 of Macaulay where 
only the projection R2 C R2 of R^ to R2 is involutive. Of course, an homogeneous system with 
constant coefficients is automatically formally integrable and one only needs to use a finite number 
of prolongations in order to obtain an involutive symbol, though it is known that 2-acyclicity is 
sufficient to obtain first order generating CC ([17]). However, it is essential to notice that it is 
only with an involutive system that we are sure that the CC system is first order both with the 
following ones in the Janet sequence. 

EXAMPLE 5.1: With = Q, m = l,n = 3, g = 2, the homogeneous second order systems 
2/33 = 0, 2/23 - 2/11 = 0, 2/22 = or 2/33 - 2/11 = 0, 2/23 = 0, 2/22 - 2/ii = both have a 2-acyclic symbol 
gs of dimension 1 at order 3 (exercise) and a trivially involutive symbol 5(4 = at order 4, such a 
result leading to only one CC of order 2 with cd{M) = 3 in both cases. We let the reader treat the 
system 2/3 = 0,2/12 = similarly and conclude (Hint: (X3,XiX2) = (X3,Xi) ^ (X3,X2) is unmixed). 
It is however not evident that the homogeneous system 2/11 = 0, 2/12 = 0, 2/13 = 0, 2/23 = of Exam- 
ple 3.33 is involutive. 

Finally, according to section 2 and 3, this first step provides the characters > ... > > 
and the smallest non-zero character a = a'^~^ ^ providing cd{M) = r, a result leading at once 
to tr{M) C M with a strict inclusion while tr-i{M) = ... = to{M) = t{M) = M. Of course, if 
a = 7^ 0, then M cannot be a torsion module and t{M) C M with a strict inclusion. The 
following example proves nevertheless that it is much more delicate to study systems with variable 
coefficients. 

EXAMPLE 5.2: With K = Q(a;2), n = 3, m = 1, g = 1, let us consider the differential module M 
defined by the trivially involutive system 2/3 — x^yi = 0. We have cd{M) = 1 but we can only say 
that cd{Dz) > l,yz £ M. If we set z = 2/2, proceeding as in Remark 2.12, we get the involutive 
system: 
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ys = x^Zi - {x^yzi 

y2 = z 

yi = zz 



2 




The differential submodule Dz C M is defined by the second order involutive system: 



■2^33 ~ 2a;^2;i3 " 


f {x-'fz.r 


= 


1 2 3 


-223 - X^Zi-2 - 


2zy 


= 


1 2 • 



and we get cd{Dz) = 1 exactly. However, even on such a very elementary example, it is not 
evident that to{M) = t{M) = M is 1-pure. We also understand that the decoupling system for 
any autonomous element in engineering sciences, like in magnetohydrodynamics, cannot be studied 
without these new techniques if wc want intrinsic results. Finally, if we denote by I the left ideal of 
D = Dy generated by 2/3 — a;^ 2/1, we notice the relation ann(G) = {xz—x^xi) — = fcid{gr{l)). 

However, we have ann{gr{Dz)) = ((xs — a;^Xi)^)X2(X3 ~ 2:^X1)) with radical equal to the prime 
ideal (xa — 2:^x1) a-s before. Hence, in this example, the strict inclusion Dz C M does not imply 
gr{{Dz) C gr{M) = G because otherwise we should get ann{G) C ann{gr{Dz) and this is the 
reason for which only the radical must be considered as it does not depend on the filtration. 



• STEP 2: Once we have obtained cd{M) = r, in order to check that M is r-pure, it remains 
to prove that tr{M) = as we already know that tr-i{M) = M. For this, the second step will 
be to use the specific properties of the Spencer form Rq+i C Ji{Rq). More generally, it is pos- 
sible to use any equivalent involutive first order system of the form i?i C Ji{E) with no zero 
order equations, that is with an induced epimorphism i?i — y E — > and such that the corre- 
sponding differential module is isomorphic and thus identified to the initial module as in Remark 
2.13 . We have now the characters a.\ > ... > a" > and the smallest non-zero character is 
is still a = a^~^ ^ providing of course the same codimension cd{M) = r as in the first step. 
Accordingly, the number r of non-zero characters and the number r of full classes is the same 
as in the previous step. However, it must be noticed that the filtration may be different and the 
following example explains once more why only the radical of the characteristic ideal must be used. 

EXAMPLE 5.3: i\: = Q, n = 2, m = 1, g = 2.For the involutive system 1/22 = 0, 2/12 = 0, the char- 
acteristic ideal is a = ((X2)^, X1X2) =^ rad{a) = (X2) =^ r = 1. Setting =y,z'^ = yi,z^ = 2/2, 
wc get the equivalent first order system ^22'^ = 0. d2z'^ = 0, d2Z^ — z^ = 0, diz^ — = 0, diz"^ = 
and the polynomial ideal generated by the 3x3 minors of the characteristic matrix is a = 
((X2)^i (X2)^Xi> X2(xi)^)- Hence the characteristic ideal is rad{a) = (X2) and r = 1 too. 

EXAMPLE 5.4: For Example 2.8 we may set z^ = y,z^ = yi,z^ = y2,z^ = 2/3 and obtain the 
first order involutive system : 



dzz'^ = 0, d^z^ = 0, d^z"^ - ^3 = 0, d^z^ -z^ = Q 12 3 
d2Z^ = Q,d2z'^ = {),d2z'^ - diz"^ = 0,d2z'^ - z^ = Q 12* 
diz"^ - z^ = 0,d-i_z^ - z'^ = Q 1 • • 

with no zero-order equation. We have a\ = 4 — 4 = 0,Q:f = 4 — 4 = 0, a} = 4 — 2 = 2 =^ r = 2 
too. We let the reader treat Example 4.2 similarly and obtain a.\ = 0, af = 0, af = 2 r = 2. 

It is at this moment that wc discover that such systems have particular properties not held by 
other systems, apart from the fact that a canonical sequence may be constructed exactly like the 
Spencer sequence or the first order part of the Janet sequence. 

Shrinking the board of multiplicative variables, we obtain from the definition of involutiveness: 

PROPOSITION 5.5: For an involutive first order system with no zero order equations and 
solved with respect to the principal {pri) first order jets expressed by means of the parametric 
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{par) other first order jets, the system obtained by looking only at the PD equations of class 1+ 
... + class i only contains di, ...,di and is still involutive Vi = 1, n, after adopting the ordering 
di+i, dn, di, di. 

EXAMPLE 5.6: Looking at Example 2.9, we notice that the systems: 



2/3 - x'^yi 


= 


4 


1 


2 3 


2/2 - 2/1 


= 


4 


1 


2 • 



{2/2-2/1=0 3 4 12 
are both involutive. Also, looking at Example 5.4, we notice that the system: 

/ d2z'^ = 0, d2Z^ = 0, d2z'^ - diz^ = 0, d2Z^ -z^ = 



diz^ 



0,diz^ 







3 


1 2 


3 


1 • 



is still involutive and we let the reader treat Example 4.2 similarly. 



We shall denote the corresponding differential module by M„_i and we notice that M — Mq 
is defined by more equations than M^-i. Accordingly, we have an epimorphism (specialization) 
Mn-i — > M — > and similarly epimorphisms Mn-i — > M„_i_i — > 0. Finally, as cd{M) = r, 
we notice that the classes n — r + 1, n arc full and we find therefore (xn)™, (Xn-r+i)™ among 
the m, X m minors with lower powers of xi, ...,Xn-r for the other minors because the numbers of 
equations of the lower classes are decreasing and thus strictly smaller than m. The characteristic 
ideal is thus (xn, Xn-r+i) if we set Xi,...,x„_r- to zero, in a coherent way. Finally, choosing 
i = n — r, we get an epimorphism — > M — y 0. The background will always indicate clearly 
the meaning of the lower index and cannot be confused with the filtration index of M. 

• STEP 3: We are now in position to study Mr with more details as a system in n — r variables 
([11], §77, p 86). Its defining system has = p = m — a equations of strict class n — r, a 

smaller number of equations of class n — r — 1, ... , and eventually an even smaller number of 
equations of class 1. Studying this system for itself, we may look for t{Mr) exactly following the 
known techniques working for any differential module, in particular double duality as described 
in section 4. However, if one could find any (relative) torsion element z e Mj., we could project 
it to an element z G M and we have N = Dz C M where we do not put a residue bar on the 
new z for simplicity. Introducing the respective annihilator ideal a and b of M and N, we should 
have C rad{a) C rad{b) as it is the only result not depending on the filtration of the modules. 
However, we know that z must be killed by at least one operator involving only di, ...,dn-r, in 
addition to the operators involving separately (dn)'" + (dn-r-i)"^ + ••• and we should obtain 

tr-i{M) = M but tr{M) ^ 0, that is M should not be pure. Hence M is r-pure if and only if Mr 
is torsion- free as a differential module over K[di, dn-r]- In such a case, the system defining Mr 
can be parametrized by a arbitrary unknowns among {y^, ...,y"^} by using a so-called minimum 
parametrization in the sense of ([24]). In actual practice, as shown on all the examples, one can use 
a relative localization with respect to (di, d„-r) only by keeping {dn-r+i, dn) untouched and 
replacing {di, dn-r) by (xi, Xn~r) considered as (constant) "parameters" in the language of 
Macaulay ([11], §43, p 45 and §77, p 86 with r instead of n — r and a different ordering). Such 
a method provides therefore a quite useful and simple test for checking purity by linking it to 
involutivity. An important intermediate result is provided by the next proposition. 

PROPOSITION 5.7: The partial localization "kill^^ the equations of class 1 up to class n — r—1 
(care) and finally only depends on the equations of strict class n — r. 

Proof. Instead of using the column n — r in the multiplicative board, we provide the proof when 
r = by using the column n. Working out as usual the first order CC, we only look at the p < m 
CC of class n for the equations of class 1 up to class n — 1 if we order the $ starting 

from the lower class involved. These p CC will be of a very specific form with a square p x p 
operator matrix for ($^, $*") with diagonal operators of the form d„ + ... where the dots denote 



21 



operators involving only in a quasi linear way with coefScicnts in K, the remaining 

of the matrix only depending on di, (i„_i for the {^'^^^ , $™) of strict class n. Therefore, if 
we have K = k = cst{K), the absolute localization is simply done by setting di — > Xi and the 
determinant of the square matrix is equal to (xnY + ••• where the dots denote a polynomial of 
degree < p in Xp with coefficients involving only xi, ■■■,Xn-i- It follows that each $p can be 

expressed as a linear combination over k{xi, Xn) of the $ of strict class n. The result is similar 
for the variable coefficient case by using the graded machinery but needs much more work. In any 
case, setting the $ of strict class n equal to zero, we should obtain a zero graded module for the 
which must be eequal to zero too. It must finally be noticed that the first order CC used 
are in reduced Spencer form as the d„ of the $ of strict class n do not appear in the CC we have 
used and these $ do not appear in the other CC too. 

Q.E.D 

EXAMPLE 5.8: Coming back to Example 2.11, we notice that the only CC is an identity in 
{u, V, w) and we may forget about u in order to obtain the new system for {v, w): 



(f)^ = d^v - x^di'v - ?; = 12 3 4 

0^ = dAW - x^diw -w = 1234 
<j)^ = dsw — x'^diw — d2V + div = 1 2 3 • 

defining a new module M with cd{M) = 1 as the class 4 is now full. The same CC as before can 
be written again in the form: 

{di - x^di - = {ds - x^dx)(j)'^ - (da - di)<?i'^ 
that we can localize in (di, ^2, ds) with: 

0^ = (^3 — x'^di)'w = {d2 — di)v V = {ds — x'^di)y, w = {d2 — d\)y 
We let the reader check that we have indeed: 

f (/)2= {di-x^di-l){d3-x'^di)y = {ds - x'^di)u = ^u = 
\ (l)^ = {di-x^di-l){d2-di)y = {d2 - di)u = =^u = 

We finally obtain for M a relative parametrization with the only constraint u = {d4—x^di — l)y = 
in a coherent way with Example 2.11. 

• STEP 4: In this section we come back to the commutative situation with a field k of constants 
and generalize the results of Macaulay in the following theorem which is recapitulating the results 

so far obtained. 

THEOREM 5.9: One has the commutative and exact diagram: 

k{Xl,-,Xn-r) ® Mr 

HXl,-,Xn-r) <8) M 
i 



Proof. For simplifying the notations in this diagram of modules over D, we have identified as a 

module over fc[di, dn-r] with Mr as a module over k[dn-r+i, c?„, di, d„-r] — k[di, dn] = 
D while the localization of M just tells that the coefficients are now in the field A;(xi, ...,Xn-r)j 
exactly following Macaulay. Moreover the central column is exact according to the definitionof Mr 
and the right column is exact because localization preserves the exactness of a sequence. 
For exampe, with fc = Q,n = 2,m = 2,(7 = 2,r = 1, the differential module M defined by 
the involutive system t/22 = 0,yi2 = may also be defined by the first order involutive system 
z'^ = y^z"^ = yi,z^ = y2 =^ d2Z^ = 0,^2^^ = 0,^22^ - z^ = 0,diz^ = 0,diz'^ ~ z"^ = 0. Then Mi 
is defined by the first order system diz^ = Q,diz^ — z^ = with torsion module t{Mi) generated 



t{Mr) 
i 

tr{M) 

i 





Mr 
i 

M 

i 
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by and the tensor product of M by k{xi) is defined by d2Z^ = 0,d2Z^ — 0,^2^^ = 0,diZ^ = 
0, z^ = 0, xi-z^ — z^ = after division by xi in a coherent way with Example 1.1. 

Q.E.D. 

Finally, when Mr is torsion- free as a differential module over k[di, ...,dn-r], then tr{M) = 
and we get the following generalization of a result provided by Macalualy ([11], §41, p 43): 

COROLLARY 5.10: The differential module M is r-pure if and only if cd{M) = r and there is 
a monomorphism — > M — > k{xi, ■■■,Xn-r) ^ M. 

• STEP 5: The final idea is to embed Mr into a free module over K[di, ...,dn-r] in order to 
parametrize the corresponding system and substitute into the equations of class n—r+1, n. How- 
ever, if we look at Example 1.2, we should find after the substitution = zi3 = 0, = Z23 = 
with one CC d2^^ — di^^ = 0, that is on one side a module L which is not 1-pure and, on the 
other side a module L having a finite free resolution with 2 operators. However, we forgot that 
M, being pure, may be identified with its embedding into its localization. Hence, we get in fact 
X1Z3 = 0, X2Z3 = and thus only 2:3 = is providing a convenient parametrizing module L. 

Our purpose is to explain and illustrate this procedure for finding such an L in the gen- 
eral situation. Again, the main idea will be provided by this example. Indeed, we obtain the 
only CC \1/ = d^^^ — ^2$^ + di^^ = 0. Substituting the paramctrization, we get of course 

= •^=^ Xiy"^ = X2y^, that is, among the two unknowns y^,y^ we are left with only one, 
say and, similarly, among the two equations we are left with only one, say because 

Xi^^ = X2^^ from the CC which is of course compatible with the localization and we choose 
2:3 = as xi-23 = Z3 = 0. 

The general situation may be treated similarly. Indeed, according to the previous step, we 
are only concerned with the equations of class n — r+1, ... , class n while the localization has 
only to do with the /3 equations of strict class n — r (care) allowing to express (3 unknowns as 
linear combinations of the a remaining unknowns with coefficients in fc(xi, Xn-r)- To each such 
equation are associated exactly r dots and each dot of index n — r + i provides a reduction of the 
respective equations of class n — r + z for z = 1, r. It follows that we are left with a equations of 
ciach such class. Whcin wc; " clclocalizc" , replacing Xi by d^, wc have to take into account the need to 
take out the denominators and may find a few "simplifications" as in the example just considered.. 
Finally, the maximum number r — 1 {care again) of dots found for one equation is obtained for the 
equations of strict class n — r+1 = n— (r — 1) and we have thus exhibited a system defining a 
module L which is r-pure and admits a free resolution with exactly (r — 1) + 1 = r operators. In 
any case, the reader must not forget that the localization of a module is useful only if we already 
know that this module is torsion-free by means of the double-duality formula t{M) = ext^{N) 
given in the introduction. 



EXAMPLE 5.11: Let M be defined by the involutive system: 

dsy'^ + diy'^ - diy^ 
dsy^ - d2y^ + diy"^ - diy^ = 
dsy^+diy"^ 
dsy'^ - diy'^ + diy^ 
d2y^ - diy^ + diy^ 
d2y^ - diy^ + diy^ 

with characters = 0, = a = 2, a} = 4. It follows that cd{M) = 1 as only the class 4 is full 
and we obtain the following relative localization showing that M is 1-pure: 

= xiy, 2/^ = xiz, y^ = (xi + X2)y, / = xiv + X2Z 

Substituting in the four equations of class 3, we only obtain the two equations: 
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r dsz + xiz =0 
1 dsy + (xi - X2)z -Xiy =0 

after a division by Xi>X2 and Xi + X2- The parametrizing module L is thus defined by the two 
equations: 

J d^z + diz =0 
\ dsy + {di - d2)z ~ diy =0 

which are differentially independent and we have the relative parametrization: 

y^=diy, y^=diz, y^ = {di+d2)y, y^ ^ diy + d2Z 

Finally, M c L ^ ass{M) c ass{L) =^ ass{M) = {{ds + di) , {ds - di)} annoiM) = 
{ds + di) D {ds — di) , a striking result showing that M can be embedded into the direct sum of two 
primary differential modules according to Remark 3.32 (See [18] for more details). 

EXAMPLE 5.12: With A; = Q, m = 1, n = 3, let us consider the polynomial map xi = u^, X2 = 

u^,X3 = "^^ in ([11], p 53). Wc have the exact sequence — > p — > k[x] — > k[u] C k{u) 
showing that p = ((X2)^(X3) - (Xi)^, (X2)^ -XiXs, (Xs)^ -X1X2) is a prime ideal ([18], p 126). The 
corresponding prime differential module M is defined by the involutive system: 



2/333 


- 2/123 


= 


1 


2 


3 


2/233 


- 2/122 


= 


1 


2 


• 


2/223 


- 2/11 


= 


1 


2 


• 


2/222 


- 2/13 


= 


1 


2 


• 


2/133 


- yii2 


= 


1 


• 


• 


2/33 - 


- 2/12 


= 


• 


• 


• 



and is 2-pure. The localized system is finite type over k{xi)[d2, ds] with par = {y, j/2, 2/3) 2/22, 2/23}- 
One can prove that is not a primary ideal even though rarf(p^) = p. 

EXAMPLE 5.13: Similarly but now with k = Q,m = l,n = 4, let us consider the polyno- 
mial map xi — uv^X2 = u,X3 = wu"^,X4 = uv^ as in ([11], p 47). We have the exact sequence 
— > p — > k[x] — > k[u,v] C k{u,v) showing that p = (x2X4 - XiX3, (xi)^ - (X2)^X3, (X3)^ - 
Xi(X4)^> (Xi)^X4 ~ X2(X3)^) is a prime ideal. It is not evident at all that the corresponding prime 
differential module M can be defined by the homogeneous involutive system (exercise): 



2/444 - 2/224 - 2/134 


- 2/123 


= 


1 


2 


3 


4 


2/344 - 2/111 




= 


1 


2 


3 


• 


2/334 - 2/114 - 2/112 




= 


1 


2 


3 


• 


2/333 - 2/124 - 2/122 


- 2/113 


= 


1 


2 


3 


• 


2/244 + 2/224 - 2/123 




= 


1 


2 


• 


• 


2/234 - 2/133 + 2/111 




= 


1 


2 


• 


• 


2/144 + 2/124 - 2/113 




= 


1 


• 


• 


• 


2/44 + 2/24 - 2/13 




= 


• 


• 


• 


• 



and is thus also 2-pure. The localized system is finite type over A:(xi, X2)[c^3; <^4] with par = 
{2/) 2/3) 2/47 2/33} and we have for example X22/34 ~ Xi2/33 + (xi)^2/ = in a coherent way with the 
comments of Macaulay in ([11], §78, p 88, formula (A) and §88,89, p 98). 

6) CONCLUSION : 

In 1916, F.S. Macaulay discovered a new localization technique for studying unmixed polyno- 
mial ideals. Wc have been able to generalize this procedure for studying pure differential modules, 
obtaining in particular a kind of relative parametrization generalizing the absolute parametriza- 
tion already known for torsion-free modules and equivalent to controllability in classical control 
theory. In the language of multidimensional systems theory, which is more intuitive, instead of 
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using arbitrary potential-like functions for the parametrization, the idea is now to use potential- 
like functions which must satisfy a kind of minimum differential constraint limiting, in some sense, 
the number of independent variables appearing in these functions, in a way similar to the situa- 
tion met in the Cartan-Khaler theorem of analysis. For such a purpose, we have exhibited new 
links between purity and involutivity, providing also a new insight into the primary decomposition 
of modules and ideals by means of tools from the formal theory of linear multidimensional systems. 

7) BIBLIOGRAPHY : 

[I] J.E. BJORK: Analytic D-modules and Applications, Kluwer, 1993. 
[2] N. BOURBAKI: Algebre homologique, chap. X, Masson, Paris, 1980. 
[3] N. BOURBAKI: Algebre commutative, chap. I-IV, Masson, Paris, 1985. 

[4] B. BUCHBERGER: Grobner bases: an Algorithmic Methods in Polynomial Ideal Theory, in: 
Recent Trends in Multidimensional System Theory, N.K. Bose ED., Reidel, Dordrecht, 1985, 184- 
232. 

[5] E. COSSERAT, F. COSSERAT: Theorie des Corps Deformables, Hermann, Paris, 1909. 

[6] W. GROBNER: Uber die Algebraischen Eigenschaften der Integrale von Linearen Differential- 

gleichungen mit Konstanten KoefRzienten, Monatsh. der Math., 47, 1939, 247-284. 

[7] M. JANET: Sur les Systemes aux derivees partielles. Journal de Math., 8, 3, 1920, 65-151. 

[8] E.R. KALMAN, Y.C. YO, K.S. NARENDA: Controllability of Linear Dynamical Systems, Con- 

trib. Diff. Equations, 1, 2, 1963, 189-213. 

[9] M. KASHIWARA: Algebraic Study of Systems of Partial Differential Equations, Memoires de 
la Societe Mathematique de France 63, 1995, (Transl. from Japanese of his 1970 Master's Thesis). 
[10] E. KUNZ: Introduction to Commutative Algebra and Algebraic Geometry, Birkhauser, 1985. 

[II] F.S. MACAULAY: The Algebraic Theory of Modular Systems, Cambridge Tracts, vol. 19, 
Cambridge University Press, London, 1916. Stechert-Hafner Service Agency, New- York, 1964. 
[12] P. MAISONOBE, C. SABBAH: D-Modules Coherents et Holonomes, Travaux en Cours, 45, 
Hermann, Paris, 1993. 

[13] D.G. NORTHCOTT: An Introduction to Homological Algebra, Cambridge University Press, 
Cambridge, 1966. 

[14] D.G. NORTHCOTT: Lessons on Rings, Modules and Muhiplicities, Cambridge University 
Press, Cambridge, 1968. 

[15] V.P. PALAMODOV: Linear Differential Operators with Constant Coefhcients, Grundlehren 
der Mathematischen Wissenschaften 168, Springer, 1970. 

[16] J.-F. POMMARET,:Dualite Differentielle et Apphcations, C. R. Acad. Sci. Paris, 320, Serie 
I, 1995, 1225-1230. 

[17] J.-F. POMMARET: Partial Differential Equations and Group Theory: New Perspectives for 
Applications, Kluwer, 1994. 

[18] J.-F. POMMARET: Partial Differential Control Theory, Kluwer, 2001. 

[19] J.-F. POMMARET: Algebraic Analysis of Control Systems Defined by Partial Differential 
Equations, in Advanced Topics in Control Systems Theory, Lecture Notes in Control and Infor- 
mation Sciences LNCIS 311, Chapter 5, Springer, 2005, 155-223. 

[20] J.-F. POMMARET: Parametrization of Cosserat Equations, Acta Mcchanica, 215, 2010, 43- 
55. 

[21] J.-F. POMMARET: Macaulay Inverse Systems Revisited, Journal of Symbolic Computation, 
46, 2011, 1049-1069. 

[22] J.-F. POMMARET: Spencer Operator and Applications: From Continuum Mechanics to 
Mathematical Physics, in "Continuum Mechanics-Progress in Fundamentals and Engineering Ap- 
phcations", Dr. Yong Gan (Ed.), ISBN: 978-953-51-0447-6, InTech, 2012, Available from: 

http://www.intechopen.com/books/continuum-mechanics-progress-in-fundamentals-and-engineerin-applications/spe: 
[23] J.-F. POMMARET: A Pedestrian Appro ach to Cosserat /Maxwell/Weyl Theory, Preprint 2012 
'http : / /hal . archives-ouvertes . fr /h al-007403"l4 
http://fr.arXiv.org/abs/1210.2675l 

[24] J.-F. POMMARET, A. QUADRAT: Localization and parametrization of linear multidimen- 
sional control systems. Systems & Control Letters 37 (1999) 247-260. 

[25] J.-F. POMMARET, A. QUADRAT: Algebraic Analysis of Linear Multidimensional Control 
Systems, IMA Journal of Mathematical Control and Informations, 16, 1999, 275-297. 



25 



[26] A. QUADRAT: Analyse Algebrique des Systemes de Controle Lineaires Multidimensionnels, 
These de Docteur de I'Ecole Nationale des Fonts et Chaussees, 1999 
( |http://www-sop. inria.fr/cafe/Alban.Quadrat/index.htnil I. 

[27] A. QUADRAT: Une Introduction a I'Analyse Algebrique Constructive et a ses Applications, 
INRIA Research Report 7354, AT-SOP Project, july 2010. Les Cours du CIRM, 1 no. 2: Journees 
Nationals de Calcul Formel (2010), p281-471 (doi:10.5802/ccirm.ll). 

[28] J.J. ROTMAN: An Introduction to Homological Algebra, Pure and Applied Mathematics, 
Academic Press, 1979. 

[29] W. M. SEILER: Involution: The Formal Theory of Differential Equations and its Applications 
to Computer Algebra, Springer, 2009, 660 pp. (See also doi:10.3842/SIGMA.2009.092 for a recent 
presentation of involution, in particular sections 3 (p 6 and reference [11], [22]) and 4). 
[30] D.C. SPENCER: Overdetermined Systems of Partial Differential Equations, Bull. Amer. 
Math. Soc, 75, 1965, 1-114. 

[31] O. ZARISKI, P. SAMUEL: Commutative Algebra, Van Nostrand, 1958. 

[32] E. ZERZ: Topics in Multidimensional Linear Systems Theory, Lecture Notes in Control and 
Information Sciences 256, Springer, 2000. 



26 



